Correlated Persistent Tunneling Currents in Glasses 
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Low temperature properties of glasses are derived within 
a generalized tunneling model, considering the motion of 
charged particles on a closed path in a double-well potential. 
The presence of a magnetic induction field B violates the time 
reversal invariance due to the Aharonov-Bohm phase, and 
leads to flux periodic energy levels. At low temperature, this 
effect is shown to be strongly enhanced by dipole-dipole and 
elastic interactions between tunneling systems and becomes 
measurable. Thus, the recently observed strong sensitivity of 
the electric permittivity to weak magnetic fields can be ex- 
plained. In addition, superimposed oscillations as function of 
the magnetic field are predicted. PACS- numbers: 61.43 Fs, 
66.35.-ha, 77.22.Ch, 03.65.Bz 

At low temperature, glasses exhibit a variety of surpris- 
ing properties of considerable theoretical interest. They 
are attributed to the existence of low-energy excitations 
present in nearly all amorphous solids and disordered 
crystals In the standard tunnehng model (STM) 
these excitations are described on a phenomenological ba- 
sis by nonintcracting two-level tunneling systems (TLS's) 
1^. Such a TLS can be approximately treated as a par- 
ticle moving in a double-well potential. At sufhciently 
low temperature only the ground states in the two wells 
are relevant, and the system is effectively restricted to 
the two-dimensional Hilbert space. Thus, the Hamilto- 
nian of a TLS has the form Hq = {\/2){/^az — Aqctj;), 
where A and Aq are the asymmetry energy and the tun- 
neling matrix element, respectively. Because of the ran- 
dom structure of glasses it is assumed that A and Aq ex- 
hibit a wide distribution according to P(A, Aq) = P/ Aq, 
where P is a constant. Treating the coupling of TLS's 
with external acoustic and electric fields as a weak per- 
turbation of Hq, the STM has successfully explained 
many of the anomalous thermal, acoustic and dielectric 
low-temperature properties of glasses. Deviations from 
the predictions of the STM |g] have been usually re- 
duced to interaction between TLS's Based on the 
spin-Boson-Hamiltonian, the elastic interaction between 
TLS's can be investigated in a nonperturbative manner 
l^. In that way it can be demonstrated that the strong 
TLS-phonon coupling essentially leads to a renormaliza- 
tion of the tunneling parameters, P and Aomin, and thus 
to the quasiparticle picture of phonon-dressed TLS's. It 
must be asked, however, whether interaction can also lead 
to an observable action of electromagnetic fields on the 
quantum-mechanical state of TLS's. This is motivated 
by the recently reported strong sensitivity of the electric 
permittivity (dielectric constant) to weak magnetic fields 



in multicomponent glasses at ultra-low temperature 
The purpose of the present paper is therefore, to investi- 
gate the influence of static electromagnetic fields on the 
energy spectrum of a TLS, and how it is modified by the 
coupling between them. 

In order to estimate electromagnetic field effects upon 
the electric permittivity of glasses in a classical approach, 
we can start from general thermodynamic principles, 
treating glasses as isotropic non-viscous dielectrics de- 
void of any free charges and free currents. The con- 
stitutive quantities of such a material can depend on 
mass density p, velocity x, temperature T, tempera- 
ture gradient VT, electric field E and magnetic flux 
density B. The general form of the constitutive equa- 
tion for a constitutive quantity such as the polarization 
P = P[p,x,T, VT, E, B] is restricted by the principle of 
material frame indifference (or objectivity, PMO) and the 
entropy principle 0. With respect to Euclidean trans- 
formations, X* = 0{t)x + b{t), where b{t) and 0{t) 
determine the translation of two frames and the rela- 
tive rotation of their axes, respectively, the PMO states 
that the constitutive equations must be invariant under 
change of frame. This means that x and E cannot occur 
separately as variables, but only in combination as elec- 
tromotive intensity £ = E + x x B. For the polarization 
the functional equation P[p, T, 0\7T, 0£, {detO)OB] = 
OP[p,T,VT,£,B] must be satisfied for all orthogonal 
matrices O. Recall that P and £ are polar vectors, while 
B is an axial vector. The general solution of the func- 
tional equation leads to the representation for the po- 
larization, whose equilibrium part (VT = 0) is given by 
P/eo = xi£ + X2{B ■ £)B + xsB x £. The susceptibil- 
ities Xa can still depend on the scalars p,T, £^ , and 
{£ ■ B)"^ . As a further restriction, it follows from the en- 
tropy principle that Xs = 0. Consequently, the electric 
permittivity e = 1 + Xi + X2BI of a rigid [x = 0) dielec- 
tric and magnetizable glass, determined by capacitance 
measurement in a plate condenser with E = {0,0, Ez), 
is expected to depend on E'^ and (and besides on 
B*, if The quadratic field dependence of e was 

indeed reported to be valid for some glasses In gen- 
eral, however, glasses are properly assumed to be lin- 
ear dielectrics, so that e(p,T) does not depend on the 
electromagnetic field. Several glasses do not show any 
measurable changes of their electric permittivity up to 
high magnetic fields and temperatures down to a few mK 
1^ . The recently observed magnetic field dependence of e 
in multicomponent glasses at ultra-low temperature 
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however, is completely different in nature from the mag- 
netoeffect in nonlinear dielectrics, and cannot be derived 
from thermodynamics assuming glasses as simple mag- 
netizable dielectrics. The knowledge of the local field 
strengths B and E only is not sufficient for the consis- 
tent description of electromagnetic field effects on the 
quantum- mechanical state of charged particles . For 
that reason the influence of a magnetic field on the energy 
spectrum of a TLS should be investigated from a micro- 
scopic point of view. In order to describe a TLS in a 
magnetic field the three-dimensional motion of a charged 
particle in an electrostatic double-well potential has to 
be considered. Imagine a hat like potential with two po- 
tential barriers in azimuthal direction along the rim of 
the hat as shown in Fig. 1. 
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FIG. 1. Hat potential. The double-well potential for a 
charged particle confined to a circular path is indicated by 
the line. The induction field B can have an arbitrary angle 
with the plane of motion 

The Hamiltonian for a non-relativistic and spinless 
charged particle of mass m and charge g, which moves 
in an uniform induction field B and in an electrostatic 
potential Y , is given by 



Ho = {p- QAf/{2m) + V. 



(1) 



The symmetric gauge for the vector potential A = 
(0,i?r/2,0) is appropriate to cylindrical co-ordinates 
(r,9,z). An uniform induction field B along the z axis 
forces the charged particle to circular motion in the per- 
pendicular plane at radii r„ with Birr^ = ipon, where 
(1)0 — h/Q is the flux quantum, and n an arbitrary inte- 
ger. The magnetic length Is = \/ h/ (QB) characterizes 
the strength of this magnetic confinement and has to be 
compared with the radius ry at which the particle is 
restrained by an electrostatic potential V{r, 9, z) of the 
form shown in Fig. 1. For weak magnetic fields, ry ^ Ib, 
the confinement force due to V exceeds the Lorentz force, 
and Landau quantization does not appear. In that case, 
the problem is reduced to that of an one-dimensional ring 
of radius r, with the Hamiltonian 



2m ^ rde 



QAof + V{e), 



(2) 



where Ag = Br/2 = <f)/{2'Kr). In the more general case 
that the induction field B has an arbitrary angle a with 
the plane of motion, the magnetic fiux through the cir- 
cular orbit of a TLS is given hy (jj = Bnr^cosa. 

For a double-well potential V{6) = Vcos{29), the 
Schrodinger equation can be reduced by means of the 
gauge transformation (fi{9) = exp{—i9(j)/(j)o)'p{d) to the 
Mathieu equation d'^(f/d9'^ + {E/a-{V/a) cos{29))(f = 0, 
where the energy E and the potential amplitude V en- 
ter in units of a = ?i^/(2mr^). Thereby, the magnetic 
flux is removed from the Hamiltonian, and the peri- 
odic boundary conditions f{9 + 27r) = y^{9) are substi- 
tuted by the magnetic flux twisted boundary conditions 
(p{9 + 2-k) — exp{—i2TT(j)/(t)Q)<f{6). Solutions with these 
boundary conditions do exist for particular values of the 
ratio E/a, which yield the energy eigenvalues. They are 
periodic in the magnetic flux (/)/(j)o as seen in Fig. 2, 
where the energies E- and E+ of the ground and the 
flrst excited state are plotted. 




FIG. 2. The two lowest energy eigenvalues as function of 
the flux ratio (j)/4>o through the circular path of a charged 
particle in a symmetric double-well potential 

For non-zero magnetic fiux the energy eigenfunctions 
composed of the odd and even Mathieu functions become 
complex. Thus, due to the magnetic induction field both 
the parity and the time reversal invariance of the TLS are 
broken. Each energy eigenstate does carry a persistent 
current 1 1 1 1 of opposite direction, given by 



I±{^)^dE±{cl>)/d<l>. 



(3) 



A net persistent tunneling current results in a magnetic 
moment of the TLS. 

For more general double- well potentials, allowing also 
an asymmetry A in the energy of the two minima, the 
harmonic approximation of the potential can be done 
around each minimum, in analogy to the STM. Then, the 
ground state and the first excited state of the TLS can 
be well approximated by a superposition of the ground 
states of each harmonic oscillator, if ^ hil > A, where 
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n is the oscillator frequency and V « (l/8)TOl7^7r^r^ is 
the potential barrier. The energy eigenvalues are found 
to be E± — {hil + A ± Eg{(j)))/2, giving the excitation 
energy (gap), 



(4) 



where t{(j)) = Aq cos(7r(/)/(/)o), with Aq w /if2exp(— A) and 
A = 2V/{hn) > 1. This is identical to the STM, when 
there the tunneling parameter Aq is substituted by the 
magnetic flux dependent tunneling splitting t(4)). 

The low-tcmpcrature thermodynamic properties cal- 
culated in the generalized model of independent TLS's, 
which are determined by the energy spectrum i?g((/)), are 
consequently periodic functions of the magnetic flux with 
a period of (/)o. Thus, the TLS energy density e((/), T) can 
be obtained by averaging the excitation energy over the 
parameters A and Aq , according to 



dAp 
Ao 



dA 



exp[Sg/(fcBT)] + l' 



(5) 



For temperatures T ^ Aomin/^B the specific heat, 
c{4>,T) = 1/p- de(4>,T) / dT , is approximately given by 



6p 



■PklT 



1 + ln 



_2fcBT_ 

imin(0) 



(6) 



where tinin(0) — Aomin ' cos(7r0/0o) is the flux depen- 
dent minimal tunneling splitting. The presence of an 
external electrical field E produces an interaction en- 
ergy of a TLS with the dipole moment p of amount 
ifint = —p ■ E. Due to E the asymmetry of the TLS 
is changed, while the change in the tunneling splitting 
can be neglected for p • E <C V. The thermodynamic 
polarization, P — —{dH/dE), can be derived from the 
effective Hamiltonian H = Hq + Hint ■ The ensemble av- 
erage over the different TLS's is done by averaging over 
possible dipole orientations and the parameters A and 
Aq with the measure P ■ dAdAo/Ag. As a result, the 
resonant part of the electric permittivity Eros becomes in 
linear response (P = eo(e — l)E), 
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dE j tanh(^). 
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^2fcBT 
(7) 



The electric permittivity depends on temperature as well 
as magnetic flux through the minimal tunneling splitting. 
There are maxima in eros('/', P) at (j)/4>o = (2n -I- l)/2, 
n — 0, f , 2, where the lower cutoff of the excitation en- 
ergy vanishes. These maxima become more pronounced 
as the temperature is lowered. Thus, the resonant part 
Ercs does depend on the magnetic field at low temperature 
where the relaxational contribution to the permittivity 
erei is negligible. Indeed, the theoretically derived depen- 
dence of the electric permittivity, Eq. (7), on magnetic 



field and temperature explains qualitatively the main fea- 
ture of the experimental results obtained in multicompo- 
nent glasses |12] as shown in Fig. 3. 
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FIG. 3. Temperature variation of the electric permittiv- 
ity Se/e = [e(T) - e(ro)]/e(To) of a BaO-Al 203-Si02 glass 
measured at a frequency of 1 kHz and a voltage of 0.75 V 
(data from Ref. (12)). 70=113 mK was taken as a reference. 
The solid circles represent the data in zero field, whereas the 
open circles represent the maximal measured increase of e 
due to a magnetic field. As shown by the solid lines, the 
data can be described by Eq.(7) in the resonant regime below 
100 mK with tmin(0 = 0)/kB = 12.2 mK (lower line) and 
imin((^ = 4'o/'2) = (upper line), respectively. Above 100 
mK the relaxational part e^d in e is additionally taken into 
account. 



In the low-temperature resonant regime, the measured 
temperature dependence of e(T, = 0) can be well de- 
scribed by Eq.(7) assuming Pp'^/eo = 1.03 • 10^^ and 
Aomin/^B = 12.2 mK. The deviation from the logarith- 
mic temperature dependence of e^cs (dielectric saturation 
|]l3| ) is lifted at magnetic fields where eres(P, = 0o/2) 
becomes maximal. In that case, the lower cutoff tmin(0) 
vanishes, and e varies logarithmically with temperature. 
The remarkable result of the experiments consists in the 
fact that both the calculated logarithmic temperature de- 
pendence of eres and its slope relative to the relaxational 
part erei of (-2):1 at higher temperature is achieved even 
at weak magnetic fields of about 0.1 T. The experimen- 
tally observed maximum in e{T,(f){B)) at P w 0.1 T, 
which slightly depends on temperature, requires to as- 
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2 • 10-1" m the TLS charge of Q « 4 • 10^|e| 



where e is the elementary charge. Before we discuss the 
origin of such a large value of Q resulting apparently 
from interactions of TLS's, we want to mention that the 
experiments described in , which were stimulated by 
our theoretical findings, have confirmed a corresponding 
flux periodic behaviour of both the specific heat and the 
electric permittivity. It should be noted that (j) in Eqs. 
(6) and (7) has to be interpreted as an effective magnetic 
flux, and an averaging over orientations and charges of 
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TLS's must be carried out in order to analyze precisely 
the oscillatory behaviour of c(T, 0) and e(T, 0) in mag- 
netic fields. 

In order to explain the large values of Q, wc have to 
consider the excitation spectrum of coupled rings. The 
Hamiltonian is given by 

i? = ^i?o. + 5]g(0»,e,)/4, (8) 

The first term is the sum of the Hamiltonians of the un- 
coupled rings, Eq. (||). The second term is the (dipole- 
dipole or elastic) interaction energy decaying with the 
distance between two TLS's i, j, and depends on 
the orientations of the dipole moments as parametrized 
by the angles 9i. We find when the interaction en- 
ergy exceeds the typical kinetic energy scale in each ring, 
that the excitation spectrum of strongly coupled rings 
equals that of one ring with an effective charge Q, which 
is the sum of the charges of the rings. This is intuitively 
clear, since then the motion of the tunneling particles 
is governed by the interaction between them, making all 
degrees of freedom massive apart from their center of 
mass motion. Thus, we have identified a possible origin 
for large values of Q. Similar situations, where the tiny 
magetic response of microscopic entities like a molecule 
is enhanced to a macroscopic magnetic field effect due to 
correlations between them, are known. One example is 
the Frederiks transition in nematic liquid crystals p5| . 
Assuming for the tunneling parameters the feasible val- 
ues P = 1045j-im-3, £;„,ax/kB = 5 K, Ao„,in/kB = IQ-^ 
K, the average distance between TLS's can be estimated 
to be r„„ — 10^^ m. Thus, a number 10^ of coupled rings 
implies mesoscopic coherence lengths, on the order of /J,m. 
For an averaged dipole moment p = 2|e| • 10""'^° m, the 
dipole-dipole interaction energy is then g/kB/r^„ ~ 100 
mK. This is exactly the temperature range in which the 
magnetic flux effects become observable. The obtained 
energy spectra for strongly coupled rings suggest the in- 
troduction of quasiparticles whose orbits are pierced by a 
flux with flux periodicity 00 = h/Q. They are excitations 
of strongly coupled TLS's with renormalized tunneling 
parameters. Due to interactions the lowest quasiparticle 
excitation energy can be considerably changed. For ex- 
ample, for two coupled rings with U > t{(j)) the splitting 
is reduced to t{(j))'^/U. 

However, there exists another phenomenon which con- 
tributes to changes in the energy spectrum and deserves 
special attention. It can be visualized by investigating in 
detail the energy levels of two interacting TLS's The 
dipole moment of an asymmetrical TLS increases with 
decreasing tunneling splitting and changes with 

variation of the external induction field. This implies 
that the dipolar coupling between two TLS's depends 
also on B, and may lead to level crossings. Energy levels 
of the coupled system may cross when their asymmetry 
is unfavourable to the interaction, and when the coupling 
increases sufficiently with magnetic field. This is the case 



if, e.g., one TLS has a small tunneling splitting and hence 
a large dipole moment. Imagine two dipole moments, 
oriented such that the dipolar interaction would flip one 
of them. When only at higher flelds the interaction is 
strong enough to achieve this, the flrst excited state and 
the ground state will cross in energy as function of B. 
Thus, a flip of a large total dipole moment of a coupled 
ring system is possible. Both level crossing and strong 
coupling result in quantum interference effects on a meso- 
scopic scale in glasses. At ultra-low temperature one may 
even expect a phase transition of coupled TLS's to occur 
as suggested in We also expect an observable action 
of the electric flux on the low-temperature dielectric re- 
sponse of glasses in alternating electrical fields [Q. In 
conclusion, we derived magnetic fiux effects in glasses by 
considering low-energy excitations as charged particles 
moving on closed paths in a double-well potential. Flux 
periodic energy levels of these tunneling states result in 
persistent tunneling currents. Due to strong coupling and 
level crossing magnetic flux effects are strongly enhanced 
below 100 mK and become measurable. 
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